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1. INTRODUCTION 

The advanced wireless communication systems promise several key performance indicators, such as 
higher data rate, low energy consumption, and higher spectral efficiency [1], which can be achieved by 
implementing beamforming capability in the transmitters and receivers. In radar applications, antennas with 
specific radiation patterns with beamforming are required [2] to ensure the radar performance to detect and 
track objects accurately. In [3], an antenna array consisting of four circularly polarized slotted elements was 
designed. Each antenna is connected to a radio frequency amplifier and a radio frequency phase shifter. In 
this way, the complex current distribution fed to each antenna can be adjusted to guarantee the correct 
beamforming. Moreover, using this approach, there are many possible errors, which degradate the quality of 
the beamforming in its sidelobe rejection and main lobe direction [4]. 

Yang et al. used low-frequency switches in all the antenna elements; in this way, they got a time- 
modulated current excitation in the antenna elements [5]-[8]. Furthermore, they used the differential 
evolution (DE) method to optimize the switching duration to get a low sidelobe by rejecting the sideband 
radiation [5]. Bio-inspired optimization methods, the ant lion optimization (ALO) in [9]-[11] and grasshopper 
optimization algorithms (GOA) in [12]-[14], were implemented to minimize the sidelobe level in a 
symmetric linear array. They delivered better results compared to well-known optimization methods such as 
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particle swarm optimization (PSO) [15], [16], firefly algorithm (FA) [17], [18], and genetic algorithm (GA) 
[19], [20]. 

The grey wolf optimization (GWO) algorithm, another bio-inspired method based on the social 
hierarchy and hunting behavior of the grey wolves, was implemented to exploit the beamforming capability 
in linear arrays [21]. Here was shown that the performance of GWO outperforms the GA. A seagull 
optimization algorithm (SOA) was recently implemented to synthesis linear arrays to obtain radiation 
patterns with low sidelobe levels with and without zeros [22]. The SOA is a new optimization method based 
on the moving and attacking behaviors of the seagull in nature. Another probabilistic optimization method 
widely used is simulated annealing, for example, designing planar multiband antennas [23] and in detecting 
of faults in linear antenna arrays [24]. 

In this paper, we would like to design multibeam arrays based on the time-modulated array [25], 
[26]. The multibeam is constructed by exploiting the resulting signals at the fundamental, first harmonic, and 
second harmonic frequencies. In contrast to other publications, we apply the algorithm of simulated 
annealing for minimizing certain cost functions. In this way we obtain the pulse duration and the switch-on 
instant in each element. Moreover, these analytical results will be verified by implementing the arrays in the 
form of microstrip arrays in a numerical simulation with HFSS. 


2. METHOD 
2.1. Time-modulated array 

Time-modulated array (TMA) is just a typical antenna array, but in each feeding line, we connect a 
switch to do a specific time variation for the complex current amplitude in each antenna [27] as given in 
Figure 1. Figure 1(a) shows a generic linear time-modulated array arranged along the z-axis. The electronic 
switches g,(t) should make the automatic beamforming possible. In order to keep the complexity of the 
practical implementation low, we restrict the time function of the switching-on and switching-off to a 
rectangular pulse function with a uniform period of T,, as given in Figure 1(b). 
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Figure 1. The principle of time-modulated array (a) TMA arrangement and (b) periodic rectangular pulse 
function of the switch 


In this condition, the array factor of the linear TMA can be seen in (1), 
AF (0d, t) = e/2nfct pee Gn(t) lh eikzncos 9 (1) 


Gn(t) is a periodic rectangular pulse function with a period T,, Tn the pulse duration, and 6,, the switch-on 
instant. This periodic switching function can be described by expanding its Fourier series as shown in (2), 


A 21T 
œ = t 
In (t) = EtL o Gnq €” “To (2) 


and especially for the periodic rectangular pulse function, the spectral contributions of all frequencies are 
given by: 
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Gag = Șn Sinc(qn§,) e Ia Gn + 20n) (3) 
where €, = Tn/To the normalized pulse duration and op = 6,,/T, the normalized switch-on instant. 


By inserting (2) and (3) into (1), we can represent the array factor of the TMA as a sum of two 
essential contributions 


AF (0, t) = [ N-i J gneikzncos | ef2nfct 


fundamental mode 
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According to (4) is the most influential formulation of TMA. The array factor is split up into the 
fundamental mode and higher-order modes, so-called harmonics, representing sideband radiation (SR). As 
we can inspect from (4), if we feed the antenna elements with the same current I, = Ip, the contribution of 
the fundamental mode is based only on the change in the new current contribution /pé,,. In this way, we can 
only have the possibility to obtain a radiation pattern with a low sidelobe level (SLL). However, if we exploit 
the harmonics, we can enhance the beamsteering capabilities of the array through the new current 
contributions with Ip ën sinc (qrën) for the amplitude and —qu(é, + 20,) for the phase. In this way, we 
have the possibility to generate additional beams in different directions as sideband radiation contributions. 


2.2. Simulated annealing 

How to adjust the switching parameters in all feeds to achieve specific radiation patterns is a kind of 
inverse problem. Therefore, solving the problems often involves optimization procedures. Although various 
methods are developed in finding of accurate and efficient meta-heuristic combinatorial optimization 
procedures, recent research hybridizes the conventional algorithms with machine learning approach [28]. 

The concept of simulated annealing originated from a physical and chemical process in metallurgy 
to obtain high-quality metals. The broad outline of the annealing process is as follows; at the beginning of the 
annealing process, a metallic structure is melted at a very high temperature, then cooled down slowly so that 
when the system reaches approximately the thermodynamic equilibrium, we have specific low-energetic 
crystalline formations. More precisely, in the beginning, at a relatively high temperature, the molecules in the 
melting metals are untidy. They are in a rapid moving condition; they move in a random way in all directions 
to achieve their appropriate states. During the cooling down process, the molecules become more ordered, 
they change their condition slowly, and the system approaches a frozen ground state at the temperature T=0. 
The process can be considered as an adiabatic thermodynamic process to the high-quality state with the 
lowest energy. Beginning with a too low initial temperature or cooling down the temperature too fast, the 
system may become quenched to form defects, or freeze out in metastable states. In metallurgy, the system is 
caught in a local minimum energy state [29]. 

The original Metropolis scheme began with a high initial temperature and a specific energy state E. 
This initial configuration is altered by keeping the temperature constant. The energy of the new proposed 
configuration may be more or lower than the previous ones. The system in metallurgy is locked in a locally 
low energy state. If the new configuration consumes less energy, the new configuration is automatically 
accepted. If the energy change is positive, the new configuration will be accepted with a probability given by 
the Boltzmann factor e~“#/". This process is repeated until the temperature is decreased and the entire 
process is repeated until a frozen condition is obtained at T=0. By analogy, the application of this Monte 
Carlo method to combinatorial problems is straighforward [30]. 

The current state of the thermodynamic system in the annealing process is analogous to the current 
solution of the combinatorial problem, the energy equation for the thermodynamic system is analogous to the 
cost function of the optimization, and the ground state is analogous to the global minimum as the problem’s 
target. The main challenge in implementing the algorithm is no obvious analogy for the temperature T 
concerning a free parameter in the combinatorial problem. Additionally, avoiding entrainment in local 
minima (quenching) is determined by how the annealing process is scheduled, i.a., the choice of initial 
temperature, the number of iterations to be conducted at each temperature, and how the temperature is 
reduced at each step as the cooling mechanism proceeds. Figure 2 depicts the flow chart of the simulated 
annealing for optimizing electromagnetic problems. 

The temperature reduction function employed in this study is T; = Ty a, with To=1000 as the initial 
temperature, T; denotes the temperature for step i, a=0.997, and the end temperature Tgyp=0.0015. 

In optimization, we define a cost function that usually must be minimized. In designing a TMA, 
some of the important goals are minimizing SLL at the fundamental mode and/or at the first harmonics, or 
minimizing the SR at the second and third harmonics, and so on (q=2, 3, ...). 
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An appropriate proposed cost function could be: 
L = wy Did, SLL y, + Wo Did, SLL pt fy + Wa Lins SRyafo ©) 
N,, Nz, and N} are the numbers of optimized SLLs/SRs at the fundamental frequencies, at the first 


harmonics, and at higher harmonics, respectively, whereby w1, W2, and w3 are the weights for these cost 
contributions. 
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Figure 2. The algorithm of simulated annealing for optimizing the desired radiation pattern 


3. RESULTS AND DISCUSSION 

In this section we design several time-modulated arrays, which fulfill some prescribed 
characteristics. Simulated annealing optimization controls the design process by varying the pulse duration 
and the switch-on instant in each feeding. The optimization goal is to reduce the cost function given in (5). 
The obtained settings are then implemented in a commercial software package (HFSS) [31]. In this way, we 
verify the target radiation diagram of each model rigorously. The antenna elements are arranged along a line 
(one-dimensional array) and are separated at a distance of 0.5. 


3.1. Optimization-based design of time-modulated arrays 
3.1.1. Time-modulated array with six elements 

In model 1, we aim to get a main beam to the direction 0 = 0° at the fundamental frequency and 
another main beam to the direction 9 = —30° at the first harmonic frequency. As model 1, we propose an 
antenna with six elements that are separated equally by a distance of A/2. With an initial start pulse duration 
and switch-on instant values, simulated annealing, according to the procedure as described in Figure 2, 
calculates the cost, which represents the deviations between the actual radiation diagram and the targeted 
ones. The optimization process continues by statistically varying the pulse duration and switch-on instant in 
each element. It ends if any stopping criterion is reached. Figure 3 summarizes the calculated result for model 
1 after the simulation process. Figure 3(a) shows the optimized pulse duration and switch-on instant of the 
switch in each feeding line. We see that the pulse durations are around 0.5, which means around 50% of the 
antennas are on. The radiation diagram, as depicted in Figure 3(b), has a main beam at 9 = 0° at the 
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fundamental frequency (red line) and, at the same time, another main beam to 0 = —30° at the first harmonic 
frequency (blue line). By inspecting the sum of the Fourier components in each antenna, the contribution of 
the first harmonic compared to the fundamental frequencies is around 1.8759/3.2948=0.57, equal t —4.9 dB. 
On the other hand, the comparison of the second harmonic to the fundamental is around 0.2991/3.2948=0.09, 
equal to —20.9 dB, which can be considered as small, so can be neglected in this case. 

In this model, we had as a target to obtain two main beams in directions V = 0° and 9 = —30°. 
However, as we can see in Figure 3(b), there are a relatively high sidelobe of the first harmonic (around —18 
dB) in the main direction of the fundamental mode (8 = 0°) and also a relatively high sidelobe of the 
fundamental mode (just around 9 dB below the main beam) in the main direction of the first harmonic (9 = 
—30°). This condition could pertube the communication link. In the following, we would like to improve this 
model by rejecting the sidelobe of the first harmonic frequency at 0 = 0° and the sidelobe of the 
fundamental frequency at 0 = —30°, while maintaining the performance of the main lobes. 

For these purposes, we modify the cost function by adding these additional targets. Figure 4 gives 
the result of the optimization process. In Figure 4(a), the switch-on duration of each antenna changes 
considerably. The central elements (#3 and #4) are switched on two times longer than the elements at the 
edges. The radiation diagram of the improved model, model 2, is given in Figure 4(b). We see that the 
sidelobe from the harmonic frequency at V = 0° reduced to more than —30 dB, and the sidelobe from the 
fundamental frequency at 0 = —30° reduced to —35 dB. In both conditions, we have more than 30 dB 
isolation to prevent any possible interference. Moreover, we observe also due to energy reduction in these 
directions, we have an enhanced energy allocation to 9 = —30° at the main beam of the first harmonic 
frequency. 


o 


4 
Element # v [°] 


(a) (b) 


Figure 3. Result for model 1 (six elements) (a) optimized pulse duration (€) and switch-on instant (0) in each 
element and (b) radiation diagram at the fundamental frequency (red), first harmonic frequency (blue) 
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Figure 4. Result for model 2 (six elements (a) optimized pulse duration (€) and switch-on instant (0) in each 
element and (b) radiation diagram at the fundamental frequency (red), first harmonic frequency (blue) 
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In this model, by inspecting the sum of the fourier components in each antenna, the contribution of 
the first harmonic compared to the fundamental frequencies is around 1.6753/2.5412=0.66, equal to —3.6 dB. 
The comparison of the second harmonic to the fundamental is around 0.176/2.5412=0.07, equal to —23.9 dB, 
which again can be considered as negligible small. 


3.1.2. Time-modulated array with sixteen elements 

In model 3, we aim to get simultaneous beamforming to 0 = 0°,9 = —45°, and 9 = 35° at the 
fundamental, first and second harmonic frequencies. At the same time, we require low signal level 
interference from each other. In order to fulfill these tight requirements, we propose an array with sixteen 
elements with equally spacing by 4/2. Figure 5 gives the pulse durations and the switch-on instants after the 
optimization process. We see, some elements are on in very short durations like elements #1, #2, #15 and #16 
as compared to others which more active in use such as the elements in the middle part of the antenna, #6, #7, 
#8, #9, #10, #11. 

The radiation diagram of these settings is given in Figure 6. We see, that the main beam at the 
fundamental frequency is oriented to 0 = 0°, where the contributions of the first and second harmonic 
frequencies are lower than —40 dB. The main beam of the first and second harmonic frequencies is directed 
to 9 = —45° and 9 = 35°, where the other spectral contributions are at around —40 dB. Overall, the possible 
interferers are lower than —30 dB than to the main beams. In this model, by inspecting the sum of the Fourier 
components in each antenna, the contribution of the first harmonic compared to the fundamental frequencies 
is around 2.9587/3.5=0.85, equal to —1.4 dB. The comparison of the second harmonic to the fundamental is 
around 1.7085/3.5=0.49, equal to —6.2 dB, as shown in Figure 6. 
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Figure 5. Optimized pulse duration (č) and switch-on instant (0) in each element of model 3 (16 elements) 
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Figure 6. Radiation diagram of model 3 (sixteen elements), fundamental frequency (red), first harmonic 
frequency (blue), second harmonic frequency (black) 
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3.2. Numerical verification 

In this section, we implement the designed array structures in microstrip technology and we 
calculate the radiation diagram numerically to verify the results obtained by the analytical approach described 
before. We use an FR4 substrate with a thickness of 1 mm, and we etch rectangular patches with a dimension 
of 11.6 mm x 11.6 cm, separated by a distance of 25 mm from each other. Figure 7 gives the computation 
model for the array with six and sixteen elements in HFSS. The antennas work at the mid-frequency of 5.8 
GHz. The patches are fed by coaxial connectors from the ground side. 


Figure 7. Models of Array with six and sixteen microstrip elements in HFSS 


Figure 8 show comparisons between the results obtained in section 3.1, illustrated as solid lines, and 
the radiation diagrams calculated by HFSS, illustrated as dashed lines. Overall, the numerical calculations 
verify the analytical results, and we see the direction of the main beams, the signal levels, and the zero 
positions are reproduced accurately both for model 1 and model 2 as given in Figure 8(a) and (b), 
respectively. Distinctive devitations are observed at the region 0 < —60° and 9 > 60°, which originate due 
to the radiation characteristics of microstrip antennas focusing the energy upwards to the region around 9 = 
0° [32]-[34]. 


Figure 8. Comparison of radiation diagram obtained by optimization (solid lines) and its verification by 
HESS (dashed lines) (a) model 1 and (b) model 2 
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Figure 9 gives the comparison for model 3. Again, the numerical results verify the obtained 
analytical design very well. We see a very small shifting of the main beam for the first harmonic (blue 
curves) and for the second harmonic (black curves) by around 2°. And we see again some deviations at the 
edges of the observed radiation region (9 < —60° and 9 > 60°) due to the radiation characteristics of the 
microstrip antennas used in numerical simulations. 


v [°] 


Figure 9. Comparison of radiation diagram obtained by optimization (solid lines) and its numerical 
verification by HFSS (dashed lines) of model 3 (sixteen elements), fundamental mode (red), first harmonic 
(blue), second harmonic (black) 


4. CONCLUSION 

In this work, we designed multibeam arrays based on the time-modulated array. The multibeam is 
contsructed by the fundamental, first harmonic, or second harmonic frequencies by varying the pulse duration 
and the switch-on instant in each element. These parameters are optimized through simulated annealing by 
minimizing the cost functions. We can show high-quality results and verify the radiation diagrams by 
numerical simulations with HFSS. Overall, the comparison of the obtained direction of the main beams, the 
signal levels, and the zero positions are very accurate. Distinctive devitations at the region 9 < —60° and 
V > 60° originate from the radiation characteristics of microstrip antennas, which focus the energy upwards 
to the region around 9 = 0°. 
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